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We review basic phenomenology on D° mixing/CP violation and recent experimental results on them. D° 
mixing is established by combining results from multiple experiments but no CP violation in the charm sector 
has been seen. D° mixing from a single experiment will clarify the size of the mixing, and observation of CP 
violation in charm decays at the present level of experimental sensitivity would be clear signal of new physics 
beyond the standard model. 



1. Introduction 



and D° states: 



Mixing of the strangeness flavor in the kaon system 
has been observed more than 50 years ago pQ and CP 
violation in the kaon system has also been well studied 
indicating, that there is no new physics (NP) beyond 
the standard model (SM) in the kaon system [2]. Also, 
the oscillation of the 6-quark flavor, in both and 
B s meson systems, has been established firmly by B- 
factories and by Tevatron and has been a leading topic 
in the flavor community of high energy physics [3]. 
CP violation effects in the B meson system have also 
been extensively carried out over last 10 years and we 
still see no strong evidence of CP violation beyond 
the SM @]. 

The study of mixing and CP violation of the D me- 
son is crucial since it involves only the up- type quark, 
which has never been studied before. Recently, mix- 
ing of D° meson has been seen by combining multiple 
experiments [5j. However, in general the intepreta- 
tion is rather difficult since the effect of the final state 
interactions are not really calculable in the SM [6]. 
Bccuase of this, the goal of search for mixing in the D 
meson system is rather to probe NP, not to constrain 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. In the 
SM, indirect CP violation in the charm system is ex- 
pected to be as small as 0(1O~ 4 ) and universal be- 
tween CP eigenstates [7J. On the other hand, the di- 
rect CP violation can be larger in SM, depending on 
the final state of the D decay of interest. In the follow- 
ing sections, we briefly discuss basic phenomenology 
of D meson mixing and CP violation and the corre- 
sponding experimental results. 



2. D° Meson Mixing 
2.1. Mixing Parameters 

In order to describe the time development of neutral 
D-meson system, one starts with writing Schrodinger 
equation for a column vector that is composed of D° 



dt 



D°(t) 



M 



(D%t)\ 

{Tf(t))> 



(1) 



where M and T are 2x2 matrices that are associ- 
ated with (D°,D ) o (D°,D ) transitions via off- 
shell (dispersive), and on-shcll (absorptive) interme- 
diate states, respectively [I]. Diagonal elements of 
the effective Hamiltonian H = M — |T are associ- 
ated with the flavor-conserving transitions, while off- 
diagonal elements are associated with flavor-changing 
transitions such as D° f-> D . The eigenstates of the 
above Schrodinger equation are parameterized as 



\Dx) oc pVl - z\D°) + qVTT^\D 
\D 2 ) oc pVl^\D°) - qVT^z\lf 



(2) 



using the notation introduced in Ref. Parameters 
p, q and z are complex- valued ones that relate flavor 
to mass eigenstates for the D-meson system. The nor- 
malized mass difference and the width difference are 
parameterized as x and y: 



IS 



(mi - m 2 )/r, 

(r x -r 2 )/2r, 



(3) 



where rrii and Ti are mass and decay rate values of 
the eigenstate \Di) (i = 1, 2). T is the average of two 
r^s. The parameters x and y are commonly called 
mixing parameters in D-meson decays and are exper- 
imentally measurable. SM calculations based on box 
diagrams alone give x ~ 10~ 5 and y ~ 10~ 7 [6], but 
are increased due to the long-distance effects. The pa- 
rameter y is dominated by long-distance effects and is 
generally considered to be insensitive to new physics. 
Therefore, x 3> y would point to NP phenomena [8]. 

In the case when D° decays to non-CP eigenstates, 
for example D° — > , one can form four different 



combinations of amplitudes as Ay 
A-j= (K-ir+\H\D°), A 



, - (K-ir+\H\D"), A f = (K-n + \H\D°), and A ? = 
{K + iv~ \H\D a ) where the first two are called "right- 
sign" decay amplitudes as they are Cabibbo favored 
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Figure 1: Projections of the proper-time distribution of 
combined D° and D wrong-sign candidates (points) and 
fit result of D° A'+tt" in (a). Fit results with and 
without mixing hypotheses are included as solid and 
dashed curves, respectively, (b) The points represent the 
difference between the data and the no-mixing fit. The 
solid curve shows the difference between fits with and 
without mixing. BaBar collaboration, Ref. |10| . 



(CF), and latter two are "wrong-sign" decay ampli- 
tudes as they are doubly Cabibbo-suppressed decays 
(DCSD) or proceed through mixing. Conventionally 
one normalizes the wrong-sign decay distributions to 
the integrated rate of right-sign decays to define r(t) 
and r(t): 



r(t) 
f(t) 



\(f\H\D°(t))\' 

\(f\n\p Q (t))f 



g+mj 1 + g-(t) 



where Xf = qAf/pAf, Xj = qAj/pAj, and g±(t) = 
i( e -izit ± e -i* 3 t) with Zl 2 _ Wl 2 /r. w 1)2 are the 
eigenvalues of Eq. . 

2.2. Semi-leptonic Decays 

Let us consider the final state of / = K + £~V(. In 
this case, Af = Aj = within the SM. The final state 
/ is only accessible through mixing and one can obtain 



r(i) 
r(f) 



\9-(t)\ 2 
\9-(t)\ 2 



-(, 2 



-(, 2 



v 2 )t 2 
y 2 )t 2 



(5) 



Note that in the SM, CP violation in charm mixing 
is small and \q/p\ ~ 1 is satisfied. Also, in the limit 
of CP conservation, r(t) = r(t). From Eq. (j5|), one 



can compute the time-integrated mixing rate relative 
to the time-integrated right-sign rate for semi-leptonic 
decays, Rm, as 



R 



1 



r(t) dt=-{x z + y z ) 



(0) 



which is a circle in x — y plane. The present world- 
average value of Rm found by the heavy flavor aver- 
aging group [5] is R M = (1-30 ± 2.69) x 1CT 4 . The 
most sensitive estimation of Rm is carried out by 
the Belle collaboration [S]. Using the decay mode of 
D*+ -> D°tt+ when D° -> K^~l + v (right-sign) or 

D° — >• D — >• K^^t^V (wrong-sign), one can infer 
the flavor of the D meson at the production by identi- 
fyig the charge of the slow pion (irf). By counting the 
yields of right-sign and wrong-sign decays one extract 
R M = (1-3 ± 2.2 ± 2.0) x 1CT 4 or Rm < 6.1 x lfr 4 at 
90% confidence level (C.L.) 0. 



2.3. Wrong-Sign Decays 

For the final state / = K + tt~ 
the ratio of decay amplitude as 




with 



one parameterizes 



O (tan 2 9 C ) (7) 



where Rd is the decay rate ratio of DCSD to CF 
modes, and 5/ is the strong phase difference between 
them. If we introduce three CP violating, real- valued 
parameters Am, Ajj, and 0, in this wrong-sign decays, 
to leading order (Ad, Am *C 1), one can write 



g+(t)Xj + g-(t) (4) r ( t ) = e - 



R D (1 + A D ) + y/R D (l + A M )(l + A D ) 
1 



x y'-t+ ^(1 + A M )R M t 2 

R D (1 - A D ) + y/R D (l-A M )(l-A D ) 



x y ' +t +-(l-A M )R M t 7 



(8) 



where y'± = y' cos 6 ± x' sin <f>, x' = x cos Sk-k + 
ysin^ifTr, and y' — y cos 8k-k — xsm Sk-w- Note that 
Sk-k is the relative strong phase between final state 
K and it, and therefore extraction of a different set 
of mixing parameters x' and y' in this final state re- 
quires the value of 5k w ■ An interference effect in the 
decay chain provides useful sensitivity to Sx-n and is 
discussed later. The BaBar collaboration looks at the 
wrong-sign decay D° —> K + n~ and fits the proper 
time distribution as shown in Fig. [ll From this, one 
extract the mixing parameters as x = (— 0.22±0.33± 
0.21) x 10~ 3 and y' = (9.7 ± 4.4 ± 3.1) x 10" 3 [TU]. 
clearly the data prefer the mixing hypothesis. 
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Figure 2: Belle collaboration 95% C.L. contours for 
(x,y): dotted (solid) corresponds to statistical (statistical 
and systematic) contour for no CP violation, and 
dash-dotted (dashed) corresponds to statistical 
(statistical and systematic) contour for the CP-allowed 
case [15| . The point is the best-fit result for no CP 
violation. 



2.4. Determination of Strong Phase 

The decay of the quantum-coherent C = — 1 state, 

-0(3770) — > D°D provides time- integrated sensitivity 
to the strong phase. The neutral charm meson in the 
CLEO-c program does not travel far enough for time- 
dependent study. Using the relations 



cos S 



K-ir 



\A{D + -> K-tt+)\ 2 - \A(D- -> K-tt+)\ 2 



where D± denotes a CP-even or -odd eigenstate, one 
can have experimental access to Sk-k- CLEO-c uses 
CP-tagged decays to obtain cosSk-k = 1-10 ± 0.35 ± 



0.07 or S K7r = 22 



,+11+ 9 
J -12-ll 



2.5. Decays to CP Eigenstates 

When D° mesons decay to CP eigenstates, for ex- 
ample / = K + K~ , there is no distinction between 
/ and /, and therefore Aj — Aj and Aj = Af. In 
this case, we define yep and A-p in order to probe the 
amount of indirect CP violation as 



Vcp 



_ " D"^K+K 



A r = 



2T 

y cos 9 — -Amx sin 9, 

T— o K+K - — ^D"^K+K- 

2T~ 



— D 



1 



Note that in the limit of CP conservation, we ex- 
pect that yep = 1 and Ay — and therefore they 
are parameters that probe CP violation phenomena 
in D°-meson decays. Substantial work on the time- 
integrated CP asymmetries in decays to CP eigen- 
states are carried out and so far all are consistent 
with no CP-violation at - 0(1)% level. The BaBar 
collaboration has studied D° — > K + K~ /ir + ir~ de- 
cays to extract a value for yep and extracted yep 
value. Experimentally, yep — (tk-tt) / (Viii) — 1, where 

(Thh) = {i~hh + Thh )/ 2 and is measured to be y C p = 
(1.16 ± 0.22 ± 0.18)%) [H] consistent with no CP vi- 
olation hypothesis. The Belle collaboration has mea- 
sured the decay-rate asymmetry for the CP-even final 
states Ap by separately determining the apparent life- 
times of D° and D in decays to the CP eigenstates 
as A T = (0.01 ± 0.30 ± 0.15)% [13,. Again, the result 
is consistent with the assumption of no CP violation 
in the charm sector. 



2.6. Dalitz Analysis 

Dalitz analysis is an invaluable technique exploited 
in many charm analyses. In general, for a three- 
body decay D — >• A + B + C, one can fully describe 
the kinematics of such decay using two parameters 
(PA + Pb) 2 and m% c = (p B + Pcf ■ They 



L AB 



AmD cos 9 — x sin 9. 



(9) 



are extremely useful since they are Lorentz invari- 
ant, the phase space of the decay is flat, and because 
of that, possible two-body resonances can be clearly 
seen. This technique has been used in light meson 
spectroscopy, CKM angle (^3(7) measurements, and 
mixing/CP violation studies. Note that this tech- 
nique can easily be extended to four and higher num- 
ber body decays. In the decay of D° —> KgTr + ir~, 
there are many quasi two-body intermediate states 
such as D° ~+ K*-n+ (CF), D Q -> K*+n- (DCSD), 
and D° — > p°Kg (CP eigenstate). Therefore, one 
form a total amplitude as 

A(m 2 _ , m 2 + ) = a r e ^ r A- (m 2 _ , m\) + a NR e^NR 

r 

A(m 2 _,ml) = J2^re^Mm-,m 2 + ) + a NR e%R 

r 

where all resonant amplitudes (A r ) are summed 
up with relative phase information as well as 
non-resonant (NR) terms. From above, time- 
dependent decay rate parameters such as e~ rt cos xTt, 
e" rt sinxrt, and e[~( 1±y ) r 'l can be extraced from the 
measurements. This is a fairly complicated analysis 
since it contains 18 quasi two-body Dalitz-plot pa- 
rameters with time-dependent unbinned maximimum 
likelihood analysis. Both Belle and BaBar collabora- 
tions analyze their data. BaBar has measured 14 x — 
(0.16 ±0.23 ±0.14)%, and y = (0.57 ± 0.20 ± 0.15)%. 
On the othe hand, Belle has measured \q/p\ and <j> 
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Figure 3: The world average contour of the mixing 
parameters x and y, indicating non-zero values of x and 

y 0- 



and Fig. [4] Even if the impact is not significant, 
the recent results from the LHCb experiment are in- 
cluded here |16| . The world average values are x = 



0.89 



+0.17 
-0.15' 



(0.63™)%, y = (0.75*0.12)%, \q/p\ 
and </>(°) = — lO.fFlg'g. These indicate that no mixing 
senario is excluded at 10ct level but there is no indica- 
tion of CP violation in the D° mixing. From this, one 
can test various NP models [17] • For example, from 
the present bound of x, one can constrain the 4th gen- 
eration quark doublet coupling as \ V u b>V c b> \ < 10 -3 for 
m b , = 500 GeV. 



3. CP Violation in D Meson Decays 

For the time-integrated search for CP violation, one 
has to extract the CP asymmetry from detector effect 
and production asymmetry in the observed asymme- 
try. To a good approximation, 



iv rec 



iv rec 



Jv rec ~ a 
— — i\, 



N D 
Jv rec 



CP 



A 



FB 



A f 



(10) 




-0.5 0.5 



x (%) 



Figure 4: The world average contour of the CP 
violation parameters in the mixing. Vertical axis is for 
and the horizontal axis is for \q/p\ [5]. 



in addition to x and y [15 



(0.80 ± 0.29; 



-0.13\ 
-0.16) 



y = (0.33±0.24t°-^)%, \q/p\ = 0.86±0.30±0.09, and 
4> = -0.24 ± 0.30 ± 0.09. Figure [2] shows the allowed 
region measured by the Belle collaboration [15] . 



2.7. World Average of Mixing Parameters 

The heavy flavor averaging group [5J average of D° 
mixing/CP violation underlying physics parameters 
from existing observables are summarized in Fig. [3] 



where N D and N D are reconstructed yields for D and 
D mesons, respectively. A^ec is the reconstructed 
asymmetry, Aqp is the CP asymmetry, Apg is the 
forward-backward asymmetry in the production, and 
A e is the one due to the charged particle reconstruc- 
tion efficiency asymmetry. To remove asymmetries 
that do not originate from CP violation, various tech- 
niques have been developed depending on availablity 
of control samples of particular decay mode of interest. 
The important point here is that to control systemat- 
ics, one really has to use real data. 

3.1. Time-integrated Search for CP 
Violation 

The CDF collaboration analyzes the data of D° — > 
K+R- and D° tt+tt- [1S : . The formulation of 
various asymmetries in this case is follows: 

A Tec (hh*) = A CP {hh) + A t {v s ) hh ' 
A ieC (Kn*) = Acp(Kit) + A e (Tr s ) K ** + A e (Kn) K "* 



\Ktt 



A V ec(KTr) = A C p(Kit) + A e (Kir) 

A C p(hh) = A Tec (hh*) - A T ec(KTT*) + A Tec (Kir) 

where (*) indicates that the D° meson flavor is tagged 
with slow pions. Note that the last equation is ob- 
tained by computing a linear combnation of above 
three. Using the relation above, Acp(tttt) = (+0.22 ± 
0.24±0.11)% and A CP {KK) = (-0.24±0.22±0.10)% 
are obtained. Both direct (a^p) and mixing-induced 
CP violation (a|?p ) contribute to the asymmetry as 



A C p = 



,ind 

l CP 



■^-aj?p where r is the mean lifetime 
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Figure 5: CP asymmetries for D + — > Kgiv + candidates 
as a function of | cos * D | in the data sample. The solid 
line represents the central value of AcP and the hatched 
region is the ±lcr interval , obtained from a minimization 
assuming no dependence on j cos 9 * D \ . 



of D° meson and (t) is the mean value of decay time 
distribution from the measurement. Note that the 
measurement errors are smaller than those from B- 
factories. 

Time-integrated CP violation search in D + — » 
Kg 7r + is carried out by both the Belle and the BaBar 
collaborations [THl [20] • In both experiments, one 
has to remove the forward-backward asymmetry from 
the production and the asymmetry from the charged 
particle reconstruction. To correct for them, Belle 
uses Df — > c/)tt + and D° — > K~ir + decay modes. 
The measured asymmetry is Aqp{D + — > A"g7r + ) = 
(-0.71 ± 0.19 ± 0.20)% where the major systemati- 
cal uncertainty is from the statistics of Df — > <j)ir + 
sample (0.18%). On the other hand, BaBar uses 
inclusive data of on- and off-resonance data for the 
correction and this enables one to reduce the sys- 
tematical uncertainty due to the correction down to 
0.08%. The corrected asymmetry value from BaBar 
is A CP (D+ -> All7r+) = (-0.44 ± 0.13 ± 0.10)%. 
The CP asymmetry from BaBar is shown in Fig. [5] 
as a function of | cos 9* D | indicating a weak deviation 
from zero. If results from two experiments are com- 
bined, one gets A C p(D+ -> K%tt+) = (-0.51±0.14)% 
and this may be the first hint of CP violation in 
the charm sector. On the other hand, this is con- 
sistent with CP violation from neutral kaon mixing 
(-0.332 ±0.006)% @]. 

A C p(D° ->• K° S P°), where P° is tt° or 77 W 



is 



measured by the Belle experiment [21] . The decay 
D* + — > D°ir+ is used in order to identify the flavor of 
the D° meson, and to correct for A n * , D° — > K~tt + 
(untagged) and D* + -> D°7r+ -> K~tt + tt+ (tagged) 



> 




E 


300- 






snts/i 




200- 


> 
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100- 
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M(hV) (GeV/c*) 



Figure 6: The invariant mass distributions of h + n 
(h + r)') in the left (right) plot after the initial selection. 
The solid histograms show -n + rj^ while the dashed 
histograms show K + r] ( -'' > final states [22] . 



are used. The results are Acp(D° —> . 
(-0.28 ± 0.19 ±0.10)%, A C p(D° -> K° s ri) = 
0.51 ± 0.16)%, and A CP {D Q -> K%i) = 
0.67 ± 0.14)%. One can assume A C p{D° 



(+0.54± 
(+0.98 ± 

K° s n°) = 

They ex- 
tract a lllu (D° K%n ) = (+0.05 ± 0.19 ± 0.10)% 
and it can be compared with the previous measure- 
ment done by the Belle experiment [13]: a' nt ^(-D° — > 
K+R-) = (-0.01±0.30±0.15)% and this is the first 



A CP 



a} n1 ^ and therefore can extract a lm ^. 
dnd^ 



experimental test of the universality of a ind in D 
decays. 

The first observations of DCSD modes D + — > 
K + rj^ is made by the Belle experiment [22]. When 
standard criteria are imposed, there is little indica- 
tion of the signal as shown in Fig. [6] In order to 
extract the such small signals from the backgrounds, 
extremely tight selection criteria are imposed based 
on a grid search technique [23] • Based on our grid 
search, a tight set of selections is imposed to the data 
and clear signals are first observed as shown in Fig. [7| 
From these, relative branching fractions are extracted 



tobe B(D+ -> K+t])/B(D 
0.14)% and B{D+ 
± 0.39 ± 0.10)%. 
which give 



n+n) = (3.06 ± 0.43 ± 
> K+i)/B{D+ -> tt+t/) = (3.77 
Using the relations in Ref. [24], 



|7T = 
|A| 2 = 

cos St a — 



3\A(K+ V )\ 2 

1 
2 

\AK + v)\ 2 
1 



\A(K+^)\ 2 + \A{K+r,')\ 2 



2\T\\A\ 
3 



2\A{K+ V )\ 2 + -\A(K+r 1 ')\ 2 



\A(K+it 



(11) 



where T (A) is the tree (annihilation) amplitude and 
A is the specified decay amplitude, and from the 
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Figure 7: The invariant mass distributions used for 
branching fraction measurements of D + — > K + w (left) 
and of D + — > K + r]' (right). Points with error bars and 
histograms correspond to the data and the fit, 
respectively |22] , 

recent branching fraction measurement of B(D + — > 
K+tt°) = (1.72 ±0.20) x 10~ 4 [25], they find that the 
relative final-state phase difference between the tree 
and annihilation in D + decays, St a, is (72 ± 9)° or 
(288 ± 9)°. This is the first experimenetal access to 
the phase difference between the tree and annihilation 
amplitudes in these decay modes. 

3.2. Time Reversal Violation 

Under the assumption of CPT invariance, T- 
violation is a signal for CP violation and therefore 
it can be studied via T-violation search. For a multi- 
particle (> 3) final state, one can form a kinematic 
product that is odd under time reveral as 

C T = pi • (p 2 x p 3 ) (12) 

where Pi is the momemtum vector of a daughter par- 
ticle i. Note that at least four different particles are 
required in the final state so that three of them are 
independent. The strong interaction dynamics can 
produce a non-zero value of the asymmetries: 

= r(c T > o) - r(c T < o) 

T ~ T(C T > 0) + T(C T < 0) 

s r(-g r >o)-r(-g r <o) 
r(-c T > o) + r(-c T < o) v ; 

where the second equation is the asymmetry for the 
CP-conjugate decay process. However, the difference 
characterizes T-violation in the weak decay process as 

At == \( a t- A T ). (14) 

This T-violation observable is measured by the BaBar 
experiment in decays of D° — > K + K~tt + tt~ and 
Df s) ->• K + KIit + tt- [26]. No evidence of T-violation 

is seen and results are At(D° —> K + K~it + tt~) = 
(1.0 ± 5.1 ± 4.4) x 10~ 3 , A T (D+ ->■ K+K%tt+it-) = 
(-12.0 ± 10.1 ± 4.6) x 10" 3 , and A T {Df -> 
K+K^n+n-) = (-13.6 ± 7.7 ± 3.4) x 10~ 3 . 



Table I Expected sensitivity with 5 ab 1 and with 50 
ab -1 in future Belle-II experiment. 



Parameters 


Present 
uncertainty 


5 ab 1 


50 ab 1 


Vcp 


±0.39 


±0.12 


±0.05 


A r 


±0.33 


±0.10 


±0.04 


x(%) 


±0.31 


±0.10 


±0.03 


y(%) 


±0.26 


±0.08 


±0.03 


\q/p\ 


±0.30 


±0.10 


±0.03 


0(rad) 


±0.30 


±0.10 


±0.03 



4. Summary and Outlook 

We reviewed present status of D° mixmg/CPV and 
D decays. The mixing, the oscillation of D° meson 
flavor, has been firmly established [5], but not from 
a single experiment yet. The search for CP violation 
in D meson decays are carried out and show no CP 
violation effect down to C(10 -3 ). There may be a hint 
of CP violation in D + — > K*gir + decay mode but it is 
consistent with CP violation from the kaon mixing. 

However, we are entering a new era of flavor physics, 
in particular for the charmed meson physics. The 
remaining yet to be analyzed data from present B- 
factories will produce even higher-sensitive results in 
near future. The LHCb is starting to demonstrate 
great ability to reconstruct charmed mesons with ex- 
tremely low background |16j . Also, planned super B- 
factories are expected to reach highest amount of data 
sample ever achieved, to 50 ab -1 [57] • Table [i] lists ex- 
pected sensitivity for various mixing and CP violation 
parameters for the Belle-II experiment, expecting to 
explore physics beyond SM, mostly through quantum 
loops in the reaction. 
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